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Abstract. 

In this contribution one examines the coupling of dark energy to the gauge fields, to 
neutrinos, and to the Higgs field. In the first case, one shows how a putative evolution 
of the fundamental couplings of strong and weak interactions via coupling to dark energy 
through a generalized Bekenstein-type model may cause deviations on the statistical nuclear 
decay Rutherford-Soddy law. Existing bounds for the weak interaction exclude any significant 
deviation. For neutrinos, a perturbative approach is developed which allows for considering 
viable varying mass neutrino models coupled to any quintessence-type field. The generalized 
Chaplygin model is considered as an example. For the coupling with the Higgs field one obtains 
an interesting cosmological solution which includes the unification of dark energy and dark 
matter. 



1. Introduction 

The nature of the dark energy and its connection with the Standard Model (SM) states is one of 
the most intriguing issues related with the negative pressure component required to understand 
the accelerated expansion of the universe. 

A putative time variation of the electromagnetic coupling on cosmological time scales from 
the observation of absorption spectra of quasars [IJ was recently reported and has led to a revival 
of interest on ideas about the variation of fundamental couplings (see Ref. [2] for a review). 
This variation is fairly natural in unification models with extra dimensions. Moreover, if there 
is a variation of the electromagnetic coupling, it is rather natural to expected variations of other 
gauge or Yukawa couplings as well. It is therefore of great interest to confirm these observations 
and to propose realistic particle physics models consistent with a cosmological variation of the 
fundamental couplings. One starts analyzing how a putative variation of the strong and weak 
couplings due to coupling to dark energy may affect the nuclear decay Rutherford-Soddy law 
[3j. In order to do that one considers the available bounds on the variation of strong and weak 
couplings and a generalized version of the Bekenstein model [1]. 
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Given that observations suggest that the variation of the electromagnetic couphng is a late 
event in the history of the universe it is natural to associate it to the observed late accelerated 
expansion of the universe [5]. This late time acceleration of the expansion can be driven by a 
cosmological constant (see [6] and references therein), by a scalar field with a suitable potential, 
quintessence [7] , or by a fluid with an exotic equation of state like the generalized Chaplygin gas 
(GCG) [Bj. For the analysis of the implications that dark energy may have on the Rutherford- 
Soddy law one considers the variation of the strong and weak couplings in the context of an 
exponential potential quintessence model and of the GCG model. 

The GCG model considers an exotic perfect fluid described by the equation of state ^ 



where A is a positive constant and a is a constant in the range < a < 1. For a = 1, the 
equation of state is reduced to the Chaplygin gas scenario. The covariant conservation of the 
energy-momentum tensor for an homogeneous and isotropic spacetime implies that 



where Ag = A/ p^i^q"^ , Pcm is the present energy density of the GCG and a the scale factor of the 
universe. Hence, one can see that at early times the energy density behaves as matter while at 
late times it behaves like a cosmological constant. This dual role is at the core of the surprising 
properties of the GCG model. Moreover, this dependence with the scale factor indicates that 
the GCG model can be interpreted as an entangled admixture of dark matter and dark energy. 

This model has been thoroughly studied from the observational point of view and it is shown 
to be compatible with the Cosmic Microwave Background Radiation (CMBR) peak location 
and amplitudes pi [10], with SNe la data [HI ttSJ [13] , gravitational lensing statistics [H], cosmic 
topology [15j, gamma-ray bursts [TB] and variation of the electromagnetic coupling [T7]. The 
issue of structure formation and its difficulties [18] have been recently addressed |19j . Most 
recent analysis based on CMBR data indicates that a < 0.25 and As > 0.93 |10j . 

One also studies the coupling of dark energy to neutrinos in the context of the so-called 
mass varying neutrino (MaVaN) models [201 [211 [23 [23 [21] • This possibility is particularly 
interesting since the coupling of neutrinos to the dark energy scalar field component may lead 
to a number of significant phenomenological consequences. If the neutrino mass nii, is generated 
by the dynamical value of a cosmologically active scalar field 4> it would be an evolving quantity. 

Hence one considers the possibility that neutrino masses arise from an interaction with the 
scalar field that drives the accelerated expansion of the universe. One considers the MaVaN 
contribution to the energy conservation equation of the cosmic fluid as a perturbation that 
iteratively modifies the background fluid equation of state [23] . 

Dark matter is most often not considered in the formulation of MaVaN models; however, 
the possibility of treating dark energy and dark matter in a unifled scheme naturally offers this 
possibility. The GCG is particularly relevant in this respect, as already mentioned, it is consistent 
with known observational constraints. Our analysis considers the coupling of MaVaN 's to the 
underlying scalar field of the GCG model. Since the neutrino contribution is perturbative, one 
obtains a small deviation from the stability condition characteristic of the unperturbed GCG 
equation of state [24j . 

Finally, one also considers a possible coupling between the Higgs boson and dark energy [25]. 
In addition to unraveling the mass generation for fermions and electroweak gauge bosons, the 
Higgs boson may also be a portal to new physics hitherto hidden in a SM singlet sector [26]. 

The most important consequences of this extension, both cosmological and phenomenological, 
have been analyzed. The phenomenological implications can arise in two ways: mixing of the 
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singlet with the Higgs boson and the possibility of invisible decay of the Higgs boson into two 
singlet bosons [26l ETJ EEl EHJ [30] . A coupling of the form ^H^H, where H is the Higgs doublet 
and <I> is the scalar singlet field, will generate a mixing between the physical Higgs boson and 
the singlet after spontaneous electroweak symmetry breaking. In addition, a quartic coupling 
like ^^H^H results in the possibility of Higgs boson decay into a pair of singlets and, in the 
case of a non-zero vacuum expectation value of the singlet, a mixing can also be induced. The 
invisible decay of the Higgs boson could be detected through the weak gauge boson fusion at 
the Large Hadron Colhder (LHC) [Slj 

From the cosmological point of view, an ultra-light scalar singlet is stable enough to be a 
good dark matter candidate, most often referred to as phion [301 [321 [331 EHl EZ] • Bounds on 
the phion-Higgs coupling can be obtained from the constraints on the phion relic abundance. 

One repeats here the analysis of Ref. [25j on the consequences of identifying the scalar singlet 
not only with the dark matter, but also with the dark energy field, responsible for the recent 
stage of accelerated expansion of the universe. It is shown that it is possible to obtain, under 
conditions, an unified picture of dark matter and dark energy, where dark energy is the zero- 
mode classical field rolling the usual quintessence potential and the dark matter candidate is the 
quantum excitation (particle) of the field, which is generated in the universe due to its coupling 
to the Higgs boson. 



2. Coupling to Gauge Fields 

As an example of the coupling of dark energy to gauge fields, one considers the impact that this 
putative coupling might have on the well established experimental fact that the nuclear decay 
rate of any nuclei is described by the statistical Rutherford-Soddy law: 

• (3) 

where X = tq being the nuclei's lifetime. For a constant A, which reflects the fact that the 
rate of transformation of an element is constant under all conditions [38], one obtains the well 
known exponential decay law: 

N{t) = iVie"^(*~*») , (4) 

where Ni is the nuclei's number counting at tj. This statistical law is a consequence of radiative 
processes of electromagnetic, strong and weak nature that take place within the nuclei. Thus, 
To is related to the amplitude of the relevant decay process being therefore a function of order 
7 of the coupling constant a^: 

To = Aa] , (5) 

with as{E) = g1{E)/hc for the strong interactions, Qs being the strong coupling constant, 
aw = GprnpC/h^ the weak coupling constant, Gp being Fermi's constant and nip the proton 
mass. Parameter A is a constant related to a specific process, suitable integration over the phase 
space, binding energy, quark masses, etc. 

On the other hand, one expects that if dark energy couples with the whole gauge sector of the 
Standard Model, this coupling can be modelled by a generalization of the so-called Bekenstein 
model ^ 

C„ = -r^f{cP)F-,F>^-^'^ , (6) 

where f{(f)) is an arbitrary function of the dark energy field, (j), and the gauge field strength. 
Given that the variation of the gauge couplings is presumably small (c.f. below), one expands 
this function to first order 
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where is the initial value of the gauge structure constant, C is a constant, (pi is the initial 
value of the quintessence field and M a characteristic mass scale of the dark energy model. We 
point out that a model with a quadratic variation is discussed in Ref. [39]. For the linear model 
it follows that 
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Thus, the gauge coupling evolution is given by 
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hence, for its variation, one obtains 
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For the electromagnetic interaction one should take into account the Equivalence Principle 
limits, which implies [4U| 



C < 7 X 10" 



(11) 



Of course, this bound does not apply to short range interactions. However, as will be seen, this 
parameter must be constrained in order to ensure that strong and weak gauge coupling do not 
change significantly so to ensure the experimental standing of the Rutherford-Soddy law. 

In what follows one considers only strong and weak decays given that the effects of a putative 
variation of the electromagnetic coupling have been the subject of various studies (see e.g. Ref. 
|41j and references therein). Thus, for strong and weak decays, keeping quark masses unchanged 
and disregarding the running of the couplings, one can write for r: 



TDE = Aficpy 



and since C << 1 
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depending on whether the process is driven by strong or weak interactions, where (po refers to 
the value of the dark energy scalar field at present. Introducing this expression into the rate of 
change of atoms 
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from which follows that [3] 



where 



N{t) ~ Ni exp [-A(At + -iCs,wI)] , 
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and At = to — ti. 

One considers now the most stringent bounds on the variation of the gauge couplings (latest 
bounds on the variation of the electromagnetic coupling can be found in [H] and references 
therein) . In what concerns the strong interaction, considerations on the stability of two-nucleon 
systems yield the bound [2] 
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On its hand, for weak interactions, considerations on the /?-decay rates of ^^'^ Re lead to the 
bound da Ha] 

^ < 3 X 10-^ . (18) 
aw 

In order to make quantitative predictions concerning the level of validity of the statistical 
Rutherford-Soddy decay law, one examines two observationally viable dark energy models. 

For simplicity, one starts considering the exponential potential V{(j)) = V^e"^*^/^^^ , in the 
case that the contribution of the scalar field is dominant, there is a family of solutions [13] 
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in terms of the scale factor of the universe. For /? < \/6, 00 can always be chosen to be zero by 
redefining the origin of 0, in which case Vq = ^(^ — l)Mp. The parameters of this model can 
be set to satisfy all phenomenological constraints (see e.g. last reference in [?])• 

One also considers the GCG model. Following Refs. |12U24j . one describes the GCG through 
a real scalar field. One starts with the Lagrangian density 



£ = -02 _ y(0) , 

where the potential for a fiat, homogeneous and isotropic universe has the following form 
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where /3 = 3(a + l)/2. 

For the energy density of the field one has 



The Friedmann equation can be written as 
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and thus 

where Q^h can be written as 
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Using Eq. ^ to obtain the pressure and integrating one obtains 
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The integral /, Eq. (I16p . has been computed for the quintessence model with an exponential 
potential (Model I) and for the GCG model (Model II) and the implications for the variation 
of the Rutherford-Soddy decay law for strong and weak interactions have been examined. For 
Model I one has set 0^ = 2 and that at present the scalar field contribution to the energy density 
is 70% of the critical density. For Model II, one has considered the values a = 0.2, Ag = 0.95, 
and used the relationship 1 + z = 1/a to express the integral / in terms of the red-shift. 

For the strong interaction decay one considers the nucleus for which tq = 4.468 x 10^ 
years. For the weak interaction decay, oone uses data from the /3~ decay of the ^^"^ Re for which 
To = 4.558 X 10^ years. Using bounds (fT7|) . and the integration interval equal to tq we 

compute C,s,w and the respective relative deviation from the Rutherford-Soddy law N/Nrs for 
the case of weak and strong interactions. Our results are shown in Table 1 for 7 = 2. 

I (yrs) \Cs\ \Cw\ I |g -1 I \wti\w -i~ 

Model I 1.22 x 10^ 7.1 x lO'^ 5.24 x lO"'^ 3.8 x IQ-^ 2.79 x IQ-'^ 
Model II 1.08 X 10^ 2.27 x IQ-^ 1.7 x lO"*^ 1.04 x 10"^ 8.1 x 10"^ 



3. Coupling to Neutrinos 

The existence of a cosmological neutrino background is a firm prediction of the cosmological 
standard model, hence any hint about this component of the universe energy density is quite 
relevant. 

The neutrino energy density and pressure are expressed through a Fermi-Dirac distribution 
function without a chemical potential, f{q), where q = and sub-index denotes present 
values. Assuming a flat Friedman- Robertson- Walker cosmology with ao = 1, one has 

P.{a,4>) = rdqq'iq' + ^I^^Y'm, (29) 

By observing that 

m,y (</>) ^q!^"^!^^^ = (Pi^ («> (A) - 3pr. (a, 0) ) , (30) 
and from Eq. (|30p . one obtains the energy- momentum conservation for the neutrino fluid 

Pu{a, (/>) + m{py{a, (j)) + py{a, </>)) = <p—^ , (31) 

where H = a/ a \s the expansion rate of the universe and the dot denotes differentiation with 
respect to cosmic time. 

The coupling between cosmological neutrinos and the scalar fleld as specified in Eq. ()30p 
is restricted to times when neutrinos are non-relativistic (NR), i. e. ^q,^""^^"^ — ni,[a) oc a^^ 
[2T| [23l SS]. On the other hand, as long as neutrinos are relativistic {Tu{a) = T^o/a » 
mi,{(j){a))), the decoupled fiuids evolve adiabatically since the strength of the coupling is 
suppressed by the relativistic increase of pressure {pi, ~ 3pu)- 



Treating the system of NR neutrinos and the scalar field as a single unified fluid (UF) which 
adiabatically expands with energy density puF = Pu + Pcf) and pressure puF = Pv + P(j> leads to 

PUF + 3H{puF + Puf) = =^ p^ + 3H{p^ + p^) = (^2) 

where the last step is derived from the substitution of Eq. ([3T]l into Eq. ([32]) . 

It is well known that the relative contribution of the energy densities components of the 
universe with respect to the one of the dark energy is on its own a problem. The assumptions 
proposed in Ref. [21] and subsequently developed elsewhere [23lll5lll6l[l8] introduce a stationary 
condition (SC) which allows circumventing the coincidence problem for cosmological neutrinos, 
by imposing that the dark energy is always diluted at the same rate as the neutrino fluid, that 
is, 

dV{(j)) dniu dpu 
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(33) 



This condition introduces a constraint on the neutrino mass since it promotes it into a dynamical 
quantity, as indicated in Eq. (I32|) . In this context, the main feature of the scenario of Ref. [21] 
is that, in what concerns to dark sector, it is equivalent to a cosmological constant-like equation 
of state and an energy density that is as a function of the neutrino mass [24j . 

At our approach, the effect of the coupling of the neutrino fluid to the scalar field fluid is 
quantified by a linear perturbation e(p (|e| << 1) such that — > (/? ~ (1 + e)(j). It then follows a 
novel equation for energy conservation 

<^ + 3ij(^ + ^KM = (34) 
dip dip dniy 

After some straightforward manipulation one obtains for the value of the coefficient of the 
perturbation [211 [26] 

drriij dpi/ 



2A. fldyw\ 

' \4> d4> ) 



for which the condition |e| << 1 is required. Upon fulfilling all known phenomenological 
requirements, the above result allows for addressing a wide class of scalar field potentials and 
related equations of state for various candidates for the dark sector (dark energy and dark 
matter), which through the SC would be incompatible with realistic neutrino mass generation 
models. 

In order to verify under which conditions Eq. ()33p agrees with the proposed perturbative 
approach for a given background equation of state, the coefficient of the linear perturbation 
should be given by 



dV{<f>) 
d4> 



d(t>\4, 



dV(<f>) \ 
d4> ) 



(|e| « 1). (36) 



This means that one must search for a neutrino mass dependence on the scale factor for which 
the above condition is satisfied. Thus, once one sets the equation of state for the dark sector, 
there will be a period at late times for which the SC and the perturbative approach match. In 
particular, this feature can be reproduced by the GCG equation of state. 

Given a potential, the explicit dependence of on the scale factor can be immediately 
obtained from Eq. ()33p . Furthermore, it is necessary to determine for which values of the scale 
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Figure 1. Present-day values of the neutrino mass tuq and tiie corresponding values of a^R 
for which the transition between the NR and UR regimes takes place in the GCG scenario with 
a = 1/2 and variable As. The increasing graylevel corresponds to increasing values of niQ, for 
the boundary values mo = 0.05 eV, 0.1 eV, 0.5 eV, leV, 5eV, respectively. 

factor the neutrino-scalar field coupling becomes important. For convenience one sets the value 
of a = a^vR for which Pu,NB. = Pu,UR holds, usually established by the condition of rriu > T^, 
that parameterizes the transition between the NR and the ultra-relativistic (UR) regime. In 
fact, this takes place when 



where x is a numerical factor estimated to be about x — 94 considering that Pu/pcru = 
mo [ey]/(94 /i^ [eT^]), where h is the value of the Hubble constant in terms of 100 km s'"^ Mpc~^ . 
Such a correspondence between qmr and mo is illustrated in the Fig. [T]for the case of a = 1/2. 
Considering the whole set of parameters that characterize the background fluid, one notices that 
it is rather difficult to see that the maximal value assumed by the e parameter corresponds to 
its present-day value. 

One observes that the interval of parameters Ag, rriQ and eventually a, for which our 
approximation can be applied (e < 1), is valid for a > unr and severally constrained by the 
condition unr ^ 1- For values of Ag (0.7 < Ag < 1) ^ [13] one finds that 0.1 < e < 1. Just 
under quite special circumstances the usual SC and the perturbative contribution of MaVaN's 
match. In the original MaVaN scenario [21] , the SC corresponds to the adiabatic solution 
(i?^ <C d^V/dc/)'^) of the scalar field equation of motion. In this case, the kinetic energy terms of 
the scalar field can be safely neglected. The consistency of our perturbative scenario with the 
stationary condition can be achieved only when the kinetic energy contribution is not relevant 
at late times. 

For As = 0, the GCG behaves always as matter, whereas for = 1, it behaves always as a 
cosmological constant. Consequently, it is natural that the relevance of the kinetic energy term 
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Figure 2. Model independent perturbative modification on squared speed of sound as a 
function of the scale factor for the neutrino-GCG coupled fluid in comparison with the adiabatic 
GCG fluid for Ag = 0.7 with a = 1, 1/2, 1/3, 10""^, for a present-day value of neutrino mass, 
mo = 0.5 eV . 

at present times is suppressed when the parameter Ag gets close to unity, which further ensures 
the agreement between the perturbative approach and the SC analysis. 

Fig. [2] illustrates the results for an increasing neutrino mass with the scale factor for a set of 
phenomenologically consistent parameters in the context of the GCG model. Interestingly, for 
mo = 0.5 ey, a fairly typical value, one can see that stable MaVaN perturbations correspond to 
a well defined effective squared speed of sound. 



The greater the mo values, the more important are the corrections to the squared speed of sound, 
up to the limit where the perturbative approach breaks down. However, one finds that as far as 
the perturbative approach is concerned, our model does not run into stability problems in the 
NR neutrino regime. In opposition, in the SC treatment, where neutrinos are just coupled to 
dark energy, cosmic expansion in combination with gravitational drag due to cold dark matter 
have a major impact on the stability of MaVaN models. Usually, for a general fluid for which the 
equation of state is known, the dominant behaviour on arises from the dark sector component 
and not by the neutrino component. For the models where the SC (cf. Eq. (I33p ) implies in a 
cosmological constant-type equation of state, = —Pfj,, one inevitably obtains = —1. 

Thus, the perturbative approach is in agreement with the assumption that the coupling 
between neutrinos and dark energy (and/or dark matter) is weak. It is found that the stability 
condition related to the squared speed of sound of the coupled fluid is predominantly governed 
by dark energy equation of state. This is rather similar to the dynamics of the weakly coupled 
cosmon fields [47J. Actually, such a troublesome behaviour should have already been observed 
as the SC implies that = — 1 from the very start, and the role of recovering the stability 
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is relegated to the neutrino contribution 1^48]. The loosening of the stationary constraint Eq. 
(I33p emerges from the dynamical dependence on 99, more concretely due to a kinetic energy 
component [21]. The knowledge of the background fluid equation of state for the dark sector 
(the GCG in our example), and the criterion for the applicability of the perturbative approach, 
do allow to overcome the negative problem, independently from the neutrino mass dependence 
set by the SC. 



4. Coupling to the Higgs field 

In order to study the coupling of dark energy to the Higgs field one considers first the coupling 
of a very light scalar to the Higgs field. One examines a model where the fields are already 
in their minimum energy configuration (contrary to the usual quintessence models, where the 
fields are displaced from their minimum). As an example one adopts the minimal spontaneously 
broken hidden sector model of Schabinger and Wells [29] involving the usual Higgs doublet H 
and a further singlet complex scalar field <I> with a potential 

V{f^,H) = -mjj\H\'^ - m||$|2 + A|i/|^ + + t]\H\'^\<^\^ . (39) 

The fields develop non- vanishing vacuum expectation values, {\H\'^) = and (I'I'P) = ^^/2, 
and in unitary gauge, the physical fields h = H — v/^/2 (the Higgs boson field) and (f) = — ^/\/2 
mix through a non-diagonal mass matrix due to a non-vanishing rj. The mixing angle lo between 
the two scalars is given by 

^""""^ " {pe - At;2) ± ^(pe^ - xv'^y^ + li^v'e ' 

In order to derive bounds on the mixing u) one considers fifth force constraints [49j. The 
mixing of the SM Higgs field with a light scalar induces a long range Yukawa-like interaction 
due to the exchange of an ultra-light singlet field with strength sin a; for electrons and qn sinu; 
for nucleons, where Qe = 2.9 x 10~^ is the electron Yukawa coupling and = 2.1 x 10~^ is the 
nucleon Yukawa coupling [50] . 

Thus, a non-relativistic test body of mass M placed in the gravitational field of the Earth 
(with mass Me) at a distance r from its center undergoes an acceleration given by: 

a = ttgr + , (41) 

where 

_ Me 

~ Mj~^ ' ^ ' 

is the usual Newtonian acceleration, Mpi is the Planck mass and 



dg%N§Nj, + gNge{N§Nl + N^N'^) + glN^Nl] (43) 



^ Mr 

( E i\ 

is the extra acceleration due to a new force arising from (j) exchange, as long as m<i.r <C 1; ^(^n'J) 
is the number of nucleons or electrons in the Earth and in the test body. 

With these definitions one finds that the difference in acceleration between two test bodies 
of distinct compositions is given by 



(44) 



where m is the average nucleoli mass and Afp is the difference in isotopic composition of the 
two test bodies: 

with N(^p^n) being the number of protons or neutrons in the test body. 

For typical isotopic differences of the order of Afp = O{10~^) and using the experimental 
result e < ©(lO^^^) one gets an estimate for upper bound on the mixing angle 

oj < 0(10-20) . (46) 

Hence the mixing of an ultra-light scalar singlet is severely constrained by fifth force-type 
experiments. However, one should point out that this is expected in the discussed model, 
since in the limit <^ v the mixing angle is given by 

\'^\^ WT- -V— (47) 
zA V ruH 

where one assumes quartic couplings of 0(1). In order for this limit to apply, the Compton 
length of the light scalar should be of the order of the Earth radius, implying in a mass of the 
order of 10"^'' eV, to be compared with a Higgs mass of the order of 100 GeV. This leads to 
10^25 _ Therefore, one concludes that there are no limits on the mixing constant r/, 
which implies that it is possible that the Higgs may have a large invisible width into very light 
particles. 

It is also possible that these light scalars are dark matter candidates. Since in this case there 
is no symmetry preventing (j) decay, one must check whether in fact they can survive till present. 
The light scalar decays through its mixing with the Higgs boson. If < 2me, the dominant 
decay rate is into photons dominated by a top quark triangle {g^it — ^9htt)^ which can be readily 
estimated: 

which implies that 

]^Q-3 

r(</> ^ 77)"' ^ —( 77^*^7 (49) 

where tu ~ lO^^s is the age of the universe. Hence, in order for the light scalar to survive till 
present one must require 

.^(^)%10- ,50) 

which is easily satisfied for the considered case. 

Hence one should consider the bounds discussed in Ref. [36], whose reasoning shows that the 
scalar (j) cannot be a dark matter candidate. In order to do that one considers the cosmological 
evolution of the light scalar field. If r] is sufficiently small, the singlet field decouples early in the 
thermal history of the universe and is diluted by subsequent entropy production. In Ref. [3l], 
one has considered out-of-equilibrium singlet production via inflaton decay in the context of 

= 1 supergravity inflationary models (see. e.g. j36j and references therein). 

On the other hand, for certain values of the coupling r/, it is possible that <I> particles are 
in thermal equilibrium with ordinary matter. In order to determine whether this is the case, 
one makes the usual comparison between the thermalization rate Tt^ and the expansion rate of 
the universe, H. In Ref. |34j it was found that if > lO^^'' the singlets can be brought into 



thermal equilibrium right after the electroweak phase transition, being therefore as abundant as 
photons. 

Since one is interested in a stable ultra-light singlet field, it is a major concern avoiding its 
overproduction if it decouples while relativistic. In fact, in this case, there is an analogue of 
Lee- Weinberg limit for neutrinos (see e.g. [53J): 

Q^h' ^ 0.08 [^^^ , (51) 

which shows that as ~ 10^^^ eV, the thermally produced ultra-light scalar (p particles cannot 
be a dark matter candidate, as for it is required that il^/i^ ~ 0.1. Of course, this model requires 
a considerable fine-tuning to generate the different scalar mass scales and it should be embedded 
in a more encompassing model, such as the Minimal Supersymmetric Standard Model with the 
addition of one singlet chiral super field [5^ . 

In what follows one considers the case where the Higgs is coupled to a quintessence field and 
pay special attention to the fact that it is not settled in its minimum. It is in this respect that 
a realistic dark energy candidate differs from the generic dark matter candidate modelled by a 
scalar <I> considered above . 

The most natural potential coupling the Higgs field to the quintessence field is of the "hybrid 
inflation" type [55] : 

y($,i7) = [/($) + A - +Ai«>2|i7p (52) 

where C/($) is the usual quintessential-type potential which can, for instance, be of the form 
of an exponential or an inverse power law in the singlet field |56[ Wl\ . The second term is the 
SM Higgs potential and the last term gives rise to Higgs-quintessence coupling, where Ai is a 
dimensionless coupling constant. 

Working in the unitary gauge and denoting h as the Higgs field, one can expand the fields 
around their classical, homogeneous vacuum configurations: 

<^{x,t) = (j)c{t) + (p{x,t); h{x,t) = -^{hc{t) +h{x,t)) (53) 

v2 

with ($(x, t)) = (pdt), {h{x, t)) = hc{t)/\/2 and the physical fields have zero vacuum expectation 
values, {if{x,t)) = {h{x,t))) = 0. 

The classical field configurations obey the equations of motion: 

dU 



+ \i<pjil = 



hc + 3Hhc + X{hi- v'')hc + 2Xi(j)ihc = 0. (54) 

Since this is a model of dynamical dark energy, the classical quintessence field has not relaxed 
to the minimum of its potential which causes the Higgs field to be displaced from the usual 
vacuum expectation value {he ^ v, where v = 246 GeV). This results in a time- varying weak 
scale or Gp, for which there are very strong bounds [58l[59] (see also the bound Eq. (fT8]l ). 

In order to guarantee that he = v, one considers the following ad hoc modifications to the last 
term of the potential Eq. (|52|) . representing the interaction between the Higgs and quintessence 
fields: 



y(cD, H)int = Xi'^Wm - v/V2f (A) 
Vi<^>,H)int = Xi{<^>-cPcf\H\^ (B) 
Vi<^,H)int = Xii<^-^ef{\H\-v/V2f {C) (55) 



For case (A) the equations of motion become: 

dU 



+ 3H(j)c + 



+ \i(l)c{hc - vY 



he + 3Hhc + \{hl - v'^)hc + V2Xi4>l{hc - v) 








(56) 



and hence he = v is a solution in which case the classical quintessence field obeys the usual 
equation: 

4>c + 3H^c + ^ 



0. 



(57) 



It is interesting to notice that this vacuum state for the $ field does not arise from a mechanism 
of spontaneous symmetry breaking, but is instead a consequence of initial conditions for the 
slow roll evolution of the quintessence field. For most of the models (pc evolves essentially as a 
logarithm of the cosmic time, so that its typical rate of change is the Hubble time: 



— H . 



(58) 



Thus, (pc evolves in a quasi-stationary regime set by the Hubble scale. 
The Lagrangian density for the physical fields is given by 



where 



C 



1 



^{d^hf - V{4>c,v,ip,h) 



(59) 



+ 



Writing 



(60) 



(61) 



and using the classical equation of motion to cancel the tadpole in Eq. (I6ip with the (pc^ term 
in Eq. ()59p one finally obtains: 



\{h^/2 + vhf - iAi(0e + ^fh^ - + 0{ip^) 



(62) 



where = U"{(l)c) and g = U"'{(l)c). 

The term in brackets in Eq. (I62p reflects the classical zero mode contribution of the scalar 
field to the energy density of the universe which causes its accelerated expansion. The other 
terms express the behaviour of the quantum excitations (particles) in the universe. In particular, 
the Higgs-quintessence interaction results in a new, time-dependent contribution to the Higgs 
boson mass given by: 



2Xv^ + Ai( 



(63) 



— — ' "-L-rc 

Since in usual quintessence models <j)c = 0{Mpi) one expects a very large deviation from 
the Standard Model value of the Higgs mass unless the coupling constant Ai is very tiny. 



Ai = 0{v'^ /Mpi) practically closing the Higgs portal. Notice, however, that unitarity in 
longitudinal gauge boson scattering is not violated since the relevant couplings are still small 
even if the Higgs boson mass is large. 

In case (B) the only difference arises in the interaction term, which becomes 



and there is no new contribution to the Higgs mass. Interestingly, in this case the coupling 
generates a new contribution to the mass of the quintessence particle 



that, contrary to usual quintessence models, does not have to be small, since it does not affect 
the evolution of 0c (i) given by Eq. ()57p . Furthermore, since the original quintessence field is 
not coupled to ordinary matter and there is no mixing to the Higgs boson, there are no bounds 
arising from fifth force constraints. Notice that in strict terms, this choice renders an ill defined 
field theory as the if mass will vary on a cosmological scale. However, for models where the 
acceleration of the expansion is transient, it is possible that U"{(pc) vanishes around present time 
(see e.g. Refs. [60l EH [62] ) . Under this condition, this case also introduces a trilinear coupling 
ip'^h which, in addition to generating an invisible Higgs boson decay, can bring the quintessence 
particle into thermal equilibrium in the early universe. The analysis of the contribution arising 
from this quintessence particle to the dark matter component follows that of Ref. [33 ^ 136 ^ 137]. An 
interesting situation occurs for the natural case of Ai = 0(0.04), which implies ~ 50 GeV. 
In this case, the particle can be brought into thermal equilibrium through the Higgs coupling 
in the early universe and, for a Higgs mass mn — 140 GeV it decouples while non-relativistic, 
generating the dark matter density in the observed range. 

Therefore, in this unified picture, dark energy is the zero-mode classical field (pdt) rolling 
down the usual quintessence potential and the dark matter candidate is the quantum excitation 
(particle) which is produced in the universe due to its coupling to the Higgs boson. 

Finally, in case (C) there are no new contributions to either to the Higgs or to the quintessence 
particle. Only a quartic interaction ip'^h'^ is generated and the quintessence particle turns out 
not to be a good dark matter candidate. 

5. Conclusions 

In this contribution, the implications of coupling dark energy to SM fields, namely gauge fields, 
neutrinos and the Higgs field, are examined. One has first considered a possible coupling of dark 
energy to gauge fields and the impact it might have on the nuclear decay Rutherford-Soddy 
law. From the available bounds on the variation of the strong and weak gauge couplings it has 
been estimated, the deviations from the Rutherford-Soddy law. It has been considered a linear 
model for the variation of the couplings. One finds that the stringent bound on the variation of 
the weak coupling arising from ^^"^Re implies in significant constraints on any variation of the 
Rutherford-Soddy law for weak interactions, actually at the level 3 x 10~^ for the quintessence 
field with an exponential potential and at 8 x 10~^ level for the GCG model. For strong 
interactions, the bounds are much less stringent and are about 4 x 10^^ for the quintessence 
model with exponential potential and of order 10~^ for the GCG model. Stronger bounds on the 
deviation of the nuclear decay law for strong interactions could be obtained if variations of the 
strong coupling were tighter. Furthermore, it is interesting that somewhat different deviations 
are found for distinct dark energy models. This may be relevant to distinguish them given that 
most often they are degenerate with respect to the observable cosmological parameters. 

In what concerns the coupling of neutrinos to dark energy, one finds that models with 
increasing neutrino mass with the scale factor is phenomenologically consistent in the context 




(64) 




(65) 



of the GCG model. For niQ = 0.5 e^, a fairly typical value, one can see that stable MaVaN 
perturbations correspond to a well defined effective squared speed of sound. The greater the 
mo values are, the more important are the corrections to the squared speed of sound, up to the 
limit where the perturbative approach breaks down [24j. 

As for the coupling of dark energy to the Higgs field, the implications of a scenario where 
the Higgs boson is coupled to a SM singlet field responsible for the accelerated expansion of the 
universe were examined. The most natural "hybrid-like" potential is disfavoured since it leads 
to a time-variation of the weak scale. A modification of the potential can however, give origin 
to the quite interesting possibility where the classical zero-mode component of the singlet field 
corresponds to the dark energy particle while its excitation plays the role of dark matter. In 
order to make this excitation consistent with the cosmological density requirement Q^h'^ ~ 0.1 
one gets a coupling to the Higgs field Ai = O(10~^), which implies = 0(10) GeV. Quite 
interestingly, this scenario can be be scrutinized in the forthcoming generation of accelerators 
through the invisible decay of the Higgs boson [25] . 
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